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After a discussion of the properties of degenerate fermion balls, we analyze 
the orbits of the star SO-1, which has the smallest projected distance to Sgr A*, 
in the supermassive black hole as well as in the fermion ball scenarios of the 
Galactic center. It is shown that both scenarios are consistent with the data, 
as measured during the last six years by Genzel et al. and Ghez et al.. We then 
consider a self-gravitating ideal fermion gas at nonzero temperature as a model 
for the Galactic halo. The Galactic halo of mass ~ 2 x lO^^M© enclosed within 
a radius of ~ 200 kpc implies the existence of a supermassive compact dark 
object at the Galactic center that is in hydrostatic and thermal equilibrium 
with the halo. The central object has a maximal mass of ~ 2.3 x lO^M© 
within a minimal radius of ~ 18 mpc or ~ 21 light-days for fermion masses 
~ 15 keV. We thus conclude that both the supermassive compact dark object 
and the halo could be made of the same weakly interacting ~ 15 keV particle. 
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1 Introduction 



In the past, self-gravitating degenerate neutrino matter has been suggested as a model 
for quasars, with neutrino masses in the 0.2keV < m < O.SMeV range M. Later it was 
used to describe dark matter in clusters of galaxies and in galactic halos, with neutrino 
masses m the 1 < m/eV < 25 range [0]. More recently, supermassive compact dark objects 
consisting of weakly interacting degenerate fermionic matter, with fermion masses in 
the 10 ^ m/keV ^ 20 range, have been proposed 0, H, ^ |^, |^ as an alternative to the 
supermassive black holes that are believed to reside at the centers of many galaxies. It has 
been pointed out that such degenerate fermion balls could cover |Q the whole range of the 
supermassive compact dark objects that have been observed so far with masses ranging 
from 10^ to 3 X 10^ Mq Most recently, it has been shown that a weakly interacting dark 
matter particle in the mass range 1 ^ m/keV ^ 5 could solve the problem of the excessive 
structure generated on subgalactic scales in A^-body and hydrodynamical simulations of 
structure formation in this Universe 0. 

So far the masses of ~ 20 supermassive compact dark objects at the center of galaxies 
have been measured using various techniques ||^. The most massive compact dark object 
ever observed is located at the center of M87 in the Virgo cluster, and it has a mass of about 
3 X IO^Mq |]T0[. If we identify this object of maximal mass with a degenerate fermion ball 
at the Oppenheimer-Volkoff (OV) limit i.e., Mqv = OMM^.^m-'^g-'^/'^ ~ 3 x IO^Mq 



0], where Mpi = ^hc/G denotes the Planck mass, this allows us to fix the fermion mass 
to m ~ 15 keV for a spin and particle-antiparticle degeneracy factor of g = 2. Such a 
relativistic object would have a radius of -Rqv = 4.45i?s — 1-5 light-days (1.19 light-days 
= 1 mpc), where -Rs is the Schwarzschild radius of the mass Moy. It would thus be 
virtually indistinguishable from a black hole of the same mass, as the closest stable orbit 
around a black hole has a radius of 3 -Rs anyway. 

Near the lower end of the observed mass range is the compact dark object located at 
the Galactic center |jl2[ with a mass of — 2.6 x IO^Mq. Interpreting this object as a 
degenerate fermion ball consisting of m ~ 15 keV and g = 2 fermions, the radius is i? ~ 21 
light-days ~ 7 x lO^Rs 0, -Rs being the Schwarzschild radius of the mass Mc. Such a 
nonrelativistic object is far from being a black hole. The shallowness of its gravitational 
potential may be probed by infalling gas and stars orbiting within this fermion ball. A 
particle at infinity at rest falling towards the fermion ball can merely acquire the escape 
velocity of ~ 1700 km/s from the center of the fermion ball. Therefore the spectrum of 
the radiation emitted by the infalling gas will be cut off at high frequencies, as is actually 
seen in the far infrared spectrum of the enigmatic radio source Sgr A* at the Galactic 
center [0. From the observed cut-off at ~ 10^'^ Hz, corresponding to emission at ~ 10 
light-days, in a thin accretion disk scenario, one can deduce a lower limit for the radius of 
-R ^ 18 light-days p. Similarly, the observed motion of stars within a projected distance 
of ~ 6 to ~ 50 light-days from Sgr A* [|12| yields, apart from the mass, an upper limit for 
the radius of the fermion ball -R ^ 22 light-days. 
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The required weakly interacting fermion of ~ 15 keV mass cannot be an active neu- 



trino, as it would overdose the Universe by orders of magnitude |T^. Moreover, an active 
neutrino of ~ 15 keV is disfavored by the experimental data on solar and atmospheric 
neutrinos, as these are most probably oscillating into active neutrinos with small ||15|| , 
and the i>e mass has been determined to be < 3 eV [l^. However, the ~15 keV fermion 
could very well be a sterile neutrino, contributing ~ 0.3 to the dark matter fraction 
of the critical density today. Indeed, as has been shown for an initial lepton asymmetry 
of ~ 10~^, a sterile neutrino of mass ~ 10 keV may be resonantly produced in the early 



Universe with near closure density, i.e. f2d ~ 1 [131 • The resulting energy spectrum of the 
sterile neutrinos is cut off for energies larger than the resonance energy, thus mimicking a 
degenerate fermion gas. As an alternative possibility, the ~ 15 keV sterile neutrino could 
be replaced by the axino or the gravitino [T^, ^ in soft supersymmetry breaking 
scenarios. 

In the recent past, galactic halos have been successfully modeled as a self-gravitating 
isothermal gas of particles of arbitrary mass, the density of which scales asymptotically 
as r~^, yielding flat rotation curves |2l|. As the supermassive compact dark objects at 
the galactic centers are well described by a gas of fermions of mass m ~ 15 keV at T = 0, 
it is tempting to explore the possibility that one could describe both the supermassive 
compact dark objects and their galactic halos in a unified way in terms of a fermion gas 
at finite temperature. We will show in this paper that this is indeed the case, and that 
the observed dark matter distribution in the Galactic halo is consistent with the existence 
of a supermassive compact dark object at the center of the Galaxy which has about the 
right mass and size, and is in thermal and hydrostatic equilibrium with the halo. 



2 Dynamics of the Stars Near the Galactic Center 

We now would like to compare the predictions of the black hole and fermion ball scenarios 
of the Galactic center, for the stars with the smallest projected distances to Sgr A*, based 
on the measurements of their positions during the last six years [0, |T^. The projected 
orbits of three stars, SO-1 (SI), SO-2 (S2) and SO-4 (S4), show deviations from uniform 
motion on a straight line during the last six years, and they thus may contain nontrivial 
information about the potential. For our analysis we have selected the star, SO-1, because 
its projected distance from Sgr A* in 1995.53, 4.4 mpc or 5.3 light-days, makes it most 
likely that it could be orbiting within a fermion ball of radius ~ 18 mpc or ~21 light-days. 
We thus may in principle distinguish between the black hole and fermion ball scenarios 
for this star. 

The dynamics of the stars in the gravitational field of the supermassive compact 
dark object can be studied solving Newton's equation of motion, taking into account the 
initial position and velocity vectors at, e.g., to = 1995.4 yr, i.e., rito) = {x,y,z) and 
r(to) = {vx, Vy, Vz). For the fermion ball the source of gravitational field is the mass Ai{r) 
enclosed within a radius r p, ^ while for the black hole it is Mc = Ai{Rc) = 2.6 x IO^Mq. 
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The X-axis is chosen in the direction opposite to the right ascension (RA), the y-axis in 
the direction of the dechnation, and the 2;-axis points towards the sun. The black hole 
and the center of the fermion ball are assumed to be at the position of Sgr A* which is 
also the origin of the coordinate system at an assumed distance of 8 kpc from the sun. 

In Figs, m and ^ the right ascension (RA) and declination of SO-1 are plotted as a 
function of time for various unobservable z's and = in 1995.4, for the black hole and 
fermion ball scenarios. The velocity components = 340 km s~^ and Vy = —1190 km 
s~^ in 1995.4 have been fixed from observations. In the case of a black hole, both RA and 
declination depend strongly on z in 1995.4, while the 2;-dependence of these quantities in 
the fermion ball scenario is rather weak. We conclude that the RA and declination data of 
SO-1 are well fitted with \z\ ~ 0.25" in the black hole scenario, and with \z\ ^ 0.1" in the 
fermion ball case (1" = 38.8 mpc = 46.2 light-days at 8 kpc). Of course, we can also try 
to fit the data varying both the unknown radial velocity and the unobservable radial 
distance z. The results are summarized in Fig. |^, where the z — v^ phase-space of 1995.4, 
that fits the data, is shown. The small range of acceptable \z\ and \vz\ values in the black 
hole scenario (solid vertical line) reflects the fact that the orbit of SO-1 depend strongly 
on z. The weak sensitivity of the orbit on z in the fermion ball case is the reason for 
the much larger z — Vz phase-space fitting the data of SO-1 |T^, as shown by the dashed 
box. The dashed and solid curves describe the just bound orbits in the fermion ball and 
black hole scenarios, respectively. The star SO-1 is unlikely to be unbound, because in 
the absence of close encounters with stars of the central cluster, SO-1 would have to fall 
in with an initial velocity that is inconsistent with the velocity dispersion of the stars at 
infinity. 

Fig. ^ shows some typical projected orbits of SO-1 in the black hole and fermion ball 
scenarios. The data of SO-1 may be fitted in both scenarios with appropriate choices of v^, 
Vy, z and Vz in 1995.4. The inclination angles of the orbit's plane, 9 = arccos (^Lz/\L\^, 

with L = mr x r, are shown next to the orbits. The minimal inclination angle that 
describes the data in the black hole case is ^ = 70°, while in the fermion ball scenario it is 
9 = 0°. In the black hole case, the minimal and maximal distances from Sgr A* are rmin 
= 0.25" and rmax = 0.77", respectively, for the orbit with z = 0.25" and Vz = which has 
a period of Tq ~ 161 yr. The orbits with z = 0.25" and Vz = 400 km s~^ or z = 0.25" 
and Vz = 700 km s~^ have periods of Tq ~ 268 yr or Tq ~ 3291 yr, respectively. In the 
fermion ball scenario, the open orbit with z = 0.1" and Vz = has a "period" of Tq ~ 
77 yr with rmin = 0.13" and rmax = 0.56". The open orbits with z = 0.1" and Vz = 400 
km s~^ or 2; = 0.1" and Vz = 900 km s~^ have "periods" of To ~ 100 yr or To ~ 1436 yr, 
respectively. 

In concluding, it is important to note that, based on the data of the star SO-1 [0, the 
fermion ball scenario cannot be ruled out. In fact, in view of the z — Vz phase space, that 
is much larger in the fermion ball scenario than in the black hole case, there is reason to 
treat the fermion ball scenario of the supermassive compact dark object at the center of 
our Galaxy with the respect it deserves. 
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3 Dark Matter in the Center and the Halo of the 
Galaxy 

Degenerate fermion balls are well understood in terms of the Thomas- Fermi theory applied 
to self-gravitating fermionic matter at T = 0. Extending this theory to nonzero 



temperature [£2], |2J], it has been shown that at some critical temperature T = 
a self-gravitating ideal fermion gas, having a mass below the OV limit enclosed in a 
spherical cavity of radius R, may undergo a first-order gravitational phase transition from 
a diffuse state to a condensed state. This is best seen plotting the energy and free energy 
as functions of the temperature which are three-valued in some temperature interval, 
exhibiting a Maxwell-Boltzmann branch at high temperatures and the degenerate branch 
at low temperatures. However, this first-order phase transition can only take place if the 
Fermi gas is able to get rid of the large latent heat which is due to the binding energy 
of the fermion ball. As the short-range interactions of the fermions are negligible, the 
gas cannot release its latent heat; it will thus be trapped for temperatures T < Tc in a 
thermodynamic quasi-stable supercooled state close to the point of gravothermal collapse. 
The Fermi gas will be caught in the supercooled state even if the total mass of the gas 
exceeds the OV limit, as a stable condensed state does not exist in this case. 

The formation of a supercooled state close to the point of gravothermal collapse, may 
be understood as a process similar to that of violent relaxation, which was introduced 
to describe rapid virialization of stars of different mass in globular clusters ||25|, ^ with- 
out invoking binary collisions of the stars, as these would not contribute significantly to 
thermalization on a scale of the age of the Universe. Through the gravitational collapse 
of a cold overdense fluctuation, ~ 1 Gyr after the Big Bang, part of gravitational energy 
transforms into the kinetic energy of random motion of small-scale density fluctuations. 
The resulting virialized cloud will thus be well approximated by a gravitationally stable 
thermalized halo. In order to estimate the mass-temperature ratio, we assume that the 
cold overdense cloud of the mass of the Galaxy M stops expanding at the time t^, reach- 
ing its maximal radius and minimal average density = 3M/(47ri?^). The total 
energy per particle is just the gravitational energy 

E^-'-^. (1) 

Assuming spherical collapse one arrives at 

Pm = —P{tm) = —^dPo{l + Zm)^, (2) 

lb lb 

where p{tm) is the background density at the time or cosmological redshift z^, and 
Po = Si^o / (8vrG) is the present critical density. We now approximate the virialized cloud 
by a singular isothermal sphere pH] of mass M and radius R, characterized by a constant 
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circular velocity B = (2T/m)^/^ and the density profile p{r) = Q"^ / A-nGr"^ . Its total energy 
per particle is the sum of gravitational and thermal energies, i.e., 

E = _i£i£ = _V. (3) 

A R A ^ ' 

Combining Eqs. ([I|), (|^), and (|^), we find 

=^-lG{m^PoM^f'\l + z^). (4) 
5 

Taking VLa = 0.3, M = 2 x IO^^Mq, = 4, and Ho = 65kms-^Mpc-\ we find 6 ~ 
220kms~^, which corresponds to the mass-temperature ratio m/T ~ 4 x 10^. 

Next, we briefly discuss the general-relativistic extension of the Thomas- Fermi theory 
23| for a self-gravitating gas of fermions with mass m and degeneracy factor g at the 



temperature T enclosed in a sphere of radius R. We denote by p, p, and n the pressure, 
energy density, and particle number density of the gas, respectively. In the following we 
use the units in which G = 1. The metric generated by the mass distribution is static, 
spherically symmetric, and asymptotically flat, i.e., 

ds'^ = ^'^dt'^ - {1 - 2M/r)-^dr^ - r\de'^ + smed(f)'^). (5) 

For numerical convenience, we introduce the parameter a = p/T and the substitution ^ = 
{ip + l)~^/^/i/m, where p is the chemical potential associated with the conserved particle 
number A^. The equation of state for a self-gravitating gas may thus be represented in 
parametric form p8| as 



1 f°° 

"""^X "^^1 + exp{[(y2 + l)V2/(<^ + 1)1/2 _ 1]^!' (6) 
^ vr^io ^l + exp{[(|/2 + l)i/2/(^ + 1)1/2 U 

v= dy 



37T^Jo ^ l + exp{[(?/2 + l)l/2/(y, + 1)1/2 _ l]c,}' W 

where appropriate length and mass scales a and b, respectively, have been chosen such 
that a = b = [2/ gY^'^ /m? . Restoring ^, c, and G, we have 
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'2 hMp\ in 2 /l5keV\' , 

^ = 1.0798 X 10^° W- km, 9 

g cm^ \ g \ m j 
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ri Ml, in /2 /l5keV\' , , 

21= 0.7251 x 10^%/- Mq. 10 

g \ g \ m 1 
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Thus fermion mass, degeneracy factor, and chemical potential are eliminated from the 
equation of state. 

Einstein's field equations for the metric @ are given by 

A 2 /ION 

Anr^p. (12) 



To these two equations we add 

dAT 



dr 



4nr'^{l -2M/r)-^^^n (13) 



dr 

imposing particle-number conservation as a condition at the boundary 

Af{R) = N. (14) 

Eqs. (|ll])-(0) should be integrated using the boundary conditions at the origin, i.e., 

ifiO) = ifo > -1 , M{0) = 0, Af{0)=0. (15) 

It is useful to introduce the degeneracy parameter r] = aip/2, which, in the Newtonian 
limit, approaches r/nr = (/^nr — V)/T, with yUnr = p — m being the nonrelativistic chemical 
potential and V the Newtonian potential. As (p is monotonously decreasing with increas- 
ing r, the strongest degeneracy is obtained at the center with rjQ = a(fo/2. The parameter 
?7o, uniquely related to the central density and pressure, will eventually be fixed by the 
requirement (0). For r > R, the function ip yields the usual empty-space Schwarzschild 
solution 

(^r =^ 1 -1, 16 

with 

M = M{R)= drAnr^pir). (17) 
Jo 

Given the temperature T, the set of self-consistency equations (|D-(|T3D, with the bound- 
ary conditions (P!^)-([T7|) defines the general-relativistic extension of the Thomas- Fermi 
equation. 



4 Numerical Results 

The numerical procedure is now straightforward. For a fixed, arbitrarily chosen a, we first 
integrate Eqs. ([TT]) and {TTj numerically on the interval [0, -R] to find the solutions for 
various central values of the degeneracy parameter rjo. Integrating (p!3D simultaneously. 
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yields J^{R) as a function of rjo. We then select the value of rjo for which Af{R) = N. The 
chemical potential /i corresponding to this particular solution is given by Eq. (|I6|) which 
in turn yields the parametric dependence on the temperature through a = fJ^/T. 

The quantities A^, T,and R are free parameters of our model and their range of values 
are dictated by the physics of the problem at hand. At T = the number of fermions 



N is restricted by the OV limit iVov = 2.89 x 10^ ^J2/ gilbkeV ImfMQ/m. However, at 
nonzero temperature, stable solutions exist with N > Nqy, depending on temperature 
and radius. In the following N is required to be of the order 2 x 10^"^ Mq/itl, so that for any 



m, the total mass is close to the estimated mass of the halo |]SD|. As we have demonstrated, 
the expected particle mass-temperature ratio of the halo is given by a ~ m/T = 4 x 10^. 
The halo radius R is in principle unlimited; in practice, however, it should not exceed 
half the average intergalactic distance. It is known that an isothermal configuration has 
no natural boundary, in contrast to the degenerate case of zero temperature, where for 
given N (up to the OV limit) the radius R is naturally fixed by the condition of vanishing 
pressure and density. At nonzero temperature, with R being unbounded, our gas would 
occupy the entire space, and fixing N would make p and p vanish everywhere. Conversely, 
if we do not fix and integrate the equations on the interval [0, oo), both M and will 
diverge at infinity for T > 0. Thus, one is forced to introduce a cutoff. In an isothermal 



model of a similar kind ||35|, the cutoff was set at the radius R, where the energy density 
was by about six orders of magnitude smaller than the central value. Our choice of 
R = 200 kpc is based on the estimated size of the Galactic halo. 

The only remaining free parameters of our model are the fermion mass m and the de- 
generacy factor g, which always appear in the combination m'^g. We fix these parameters 
at m = 15 keV and g = 2, and justify this choice a posteriori. 

The results of numerical integration of Eqs. ( pH]) and (|T^), without restricting A^, are 
presented in Fig. ^, where we plot the particle number A^ as a function of the central 
degeneracy parameter r]o for several values of a close to 4 x 10^. For fixed A^, there is a 
range of a, where the Thomas-Fermi equation has multiple solutions. For example, for 
A^ = 2 X 10^^ and a = 4 x 10^ six solutions are found, which are denoted by (1), (2), 
(3), (3'), (2'), and (!') corresponding to the values r/o = -30.528, -25.354, -22.390, 29.284, 
33.380, and 40.479, respectively. In Figs. || and ^ we plot the corresponding density 
profiles and enclosed masses, respectively. For negative central value 770, for which the 
degeneracy parameter is negative everywhere, the system behaves basically as a Maxwell- 
Boltzmann isothermal sphere. Positive values of the central degeneracy parameter r]o are 
characterized by a pronounced central core of mass of about 2.5 x IO^Mq within a radius 
of about 20 mpc. The presence of this core is obviously due to the degeneracy pressure of 
the Fermi-Dirac statistics. A similar structure was obtained in coUisionless stellar systems 



modeled as a nonrelativistic Fermi gas 



Figs. H and show two important features. First, a galactic halo at a given temper- 
ature T may or may not have a central core depending whether the central degeneracy 
parameter rjo is positive or negative. Second, the closer to zero t]q is, the smaller the 
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radius is at which the asymptotic behavior of the density begins. The flattening of 
the Galactic rotation curve begins in the range 1 ^ r/kpc ^ 10, hence the solution (3') 
most likely describes the Galaxy's halo. This may be verified by calculating the rotational 
curves in our model. We know already from the estimate that our model yields the 
correct asymptotic circular velocity of 220 km/s. In order to make a more realistic com- 
parison with the observed Galactic rotation curve, we must include two additional matter 
components: the bulge and the disk. The bulge is modeled as a spherically symmetric 
matter distribution of the form [ plU 



g-/is „oo ^—hsu 



2s3 Jo [(m+1)8-1]V2' ^ ^ 

where s = (r/ro)^/4, tq is the effective radius of the bulge and h is a parameter. We 
adopt To = 2.67 kpc and h yielding the bulge mass Mb = 1.5 x IO^^Mq In Fig. ^ the 



mass of halo and bulge enclosed within a given radius is plotted for various tiq. Here, the 
gravitational backreaction of the bulge on the fermionic halo has been taken into account. 
The data points, indicated by squares, are the mass Mc = 2.6 X IO^Mq within 18 mpc, 
estimated from the motion of the stars near Sgr A* [jl2l, and the mass M50 = 5.4l3;g x 10^^ 
within 50 kpc, estimated from the motions of satellite galaxies and globular clusters . 
Variation of the central degeneracy parameter r/o between 24 and 32 does not change the 
essential halo features. 

In Fig. 1^ we plot the circular velocity components of the halo, the bulge, and the disk. 
The contribution of the disk is modeled as 

where we take Tq = 13.5 kpc and Od = 100 km/s. Here it is assumed for simplicity that 
the disk does not influence the mass distribution of the bulge and the halo. Choosing 
the central degeneracy rjo = 28 for the halo, the data by Merrifield and Oiling [|3^ are 
reasonably well fitted. 

We now turn to the discussion of our choice of the fermion mass m = 15 keV for a 
degeneracy factor g = 2. To this end we investigate how the mass of the central object, 
i.e., the mass Mc within 18 mpc, depends on m in the interval 5 to 25 keV, for various tjq. 
We find that m ^ 15 keV gives always the maximal value of ranging between 1.7 to 
2.3 xIO^Mq for 7]q between 20 and 28. Hence, with m ~ 15 keV we get the value closest 
to the mass of the central object estimated from the motion of the stars near Sgr A* 

We now present the results of the calculations for fixed particle number and temper- 



atures near the point of gravothermal collapse. In Fig. |TU| the energy per particle defined 
as E = M/N — m is plotted as a function of temperature for fixed N = 2 x 10^^. The 



plot looks very much like that of a canonical Maxwell-Boltzmann ensemble [^, with 
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one important difference: in the Maxwell-Boltzmann case the curve continues to spiral 
inwards ad infinitum approaching the point of the singular isothermal sphere, that is 
characterized by an infinite central density. In Fermi-Dirac case the spiral consists of 
two almost identical curves. The inwards winding of the spiral begins for some negative 
central degeneracy and stops at the point T = 2.3923 x 10~''m, E = —1.1964 x 10~'^m, 
where 7]q becomes zero. This part of the curve, which basically depicts the behavior of 
a nondegenerate gas, we call Maxwell-Boltzmann branch. By increasing the central de- 
generacy parameter further to positive values, the spiral begins to unwind outwards very 
close to the inwards winding curve. The outwards winding curve will eventually depart 
from the Maxwell-Boltzmann branch for temperatures T ^ 10~^m. Further increase of 
the central degeneracy parameter brings us to a region, where general-relativistic effects 
become important. The curve will exhibit another spiral for temperatures and energies 
of the order of a few 10~^m approaching the limiting temperature Too = 2.4 x 10~^m and 
energy E^q = 3.6 x 10~^m with both the central degeneracy parameter and the central 
density approaching infinite values. It is remarkable that gravitationally stable configura- 
tions with arbitrary large central degeneracy parameters exist at finite temperature even 
though the total mass exceeds the OV limit by several orders of magnitude. 



5 Conclusions 

In summary, using the Thomas-Fermi theory, we have shown that a weakly interacting 
fermionic gas at finite temperature yields a mass distribution that successfully describes 
both the center and the halo of the Galaxy. For a fermion mass m ~ 15 keV, a reasonable 
fit to the rotation curve is achieved with the temperature T = 3.75 meV and the degen- 
eracy parameter at the center rjQ = 28. With the same parameters, we obtain the mass 
M50 = 5.04 X W^Mq and M200 = 2.04 x W'^Mq within 50 and 200 kpc, respectively. 
These values agree quite well with the mass estimates based on the motions of satellite 
galaxies and globular clusters |^ . Moreover, the mass of ~ 2.27 x IO^Mq, enclosed 



within 18 mpc, agrees reasonably well with the observations of the compact dark object 
at the center of the Galaxy. We thus conclude that both the Galactic halo and center 
could be made of the same fermions. 

An observational consequence of this unified scenario of fermion ball and fermion halo 
at finite temperature could be the direct observation of the radiative decay of the fermion 
(assumed here to be a sterile neutrino) into a standard neutrino, i.e., / w-f. The X-ray 
luminosity of the compact dark object is most easily observed. If the lifetime for the decay 
/ — ^ z/7 is 0.82 xlO^^ yr, the luminosity of a = 2.6 x IO^Mq fermion ball would be 0.9 
xlO'^^ erg s~^. This is consistent with the upper limit of the X-ray luminosity of ~ (0.5 
- 0.9) xlO^^ erg s~^ of the source with radius 0.5" ~ 23 light-days, whose center nearly 
coincides with Sgr A*, as seen by the Chandra satellite in the 2 to 7 keV band The 
lifetime is proportional to sin~^ 9, 9 being the unknown mixing angle of the sterile with 
active neutrinos. With a lifetime of 0.82 x 10^^ yr we obtain an acceptable value for the 
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mixing angle squared of = 1.4 x 10^^^. The X-rays originating from such a radiative 
decay would contribute about two orders of magnitude less than the observed diffuse X- 
ray background at this wavelength if the sterile neutrino is the dark matter particle of 
the Universe. The signal observed at the Galactic center would be a sharp X-ray line at 
~ 7.5 keV for g = 2 and ~ 6.3 keV for (7 = 4. This line could be misinterpreted as the 
Fe Ka line at 6.67 keV. Scattering with baryonic matter within the Galactic center could 
distribute the energy more evenly in the 2 to 7 keV band. The X-ray luminosity would 
be tracing the fermion matter distribution, and it could thus be an important test of the 
fermion ball scenario. Of course the angular resolution would need to be ^ 0.1" and the 
sensitivity would have to extend beyond 7 keV. 
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Figure 1: Right ascension of SO-1 versus time for various \z\ and Vx = 340 km s ^, 



Vy = -1190 km s~^ and i;^ = in 1995.4. 



Figure 2: Declination of SO-1 versus time for various \z\ and Vx = 340 km s \ Vy = —1190 
km s~^ and = in 1995.4. 



Figure 3: The z — phase-space that fits the SO-1 data. 
Figure 4: Examples of typical orbits of SO-1. 



Figure 5: Number of particles versus central degeneracy parameter for m/T = 4 x 10^ 
(solid), 3.5 X 10^ (short dashs), 4.5 x 10^ (long dashs), and 5 x 10^ (dot-dashed line). 



Figure 6: The density profile of the halo for a central degeneracy parameter 770 = 
(dotted line) and for the six r/o- values discussed in the text. Configurations with negative 
r/o ((l)-(3)) are depicted by the dashed and those with positive r/o ((l')"(3')) by the solid 
line. 



Figure 7: Mass of the halo Mh(r) enclosed within a radius r for various central degeneracy 
parameters tiq as in Fig. 0. 



Figure 8: Enclosed mass of halo plus bulge versus radius for rjo = 24 (dashed), 28 (solid), 
and 32 (dot-dashed line). 



Figure 9: Fit to the rotation curve of the Galaxy. The data points are from for 
Ro = 8.5 kpc and Go = 220 km/s. 



Figure 10: Energy (shifted by 12 x 10 ^m) versus temperature (shifted by —24 x 10 ^m) 
both in units of 10-^°m, for fixed iV = 2 x W^Mq/ui 
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